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Abstract 

The cross section for inclusive charm production by a low-energy antiproton beam is calculated 
using the parton model and next-to-leading order perturbative QCD. For an antiproton beam with 
a momentum of 15 GeV, the charm cross section at next-to-leading order in the QCD coupling 
constant changes by more than an order of magnitude as the charm quark mass is varied from 1.3 
to 1.7 GeV. The variations can be reduced by demanding that the same value of the charm quark 
mass give the measured charm cross sections for fixed-target experiments with a proton beam. The 
resulting estimate for the charm cross section from a low-energy antiproton beam is large enough 
to allow the study of charm meson mixing. 

PACS numbers: 12.38.-t, 12.39. St, 13.20. Gd, 14.40. Gx 



Flavor physics in the charm sector is an aspect of elementary particle physics that has 
not been fully explored. The Standard Model of particle physics predicts the rates for charm 
meson mixing and for CP violation in charm meson decays to be very small, which makes 
the charm sector a promising window into new physics beyond the Standard Model [lj]. 
Mixing of the neutral charm mesons D° and D° was finally observed by several experiments 
in 2007 [1, 0, [5[ . Searches for CP violation in the charm sector would require significantly 
larger data samples. One possible source of large data samples of flavor-tagged charm mesons 
is a low-energy antiproton beam, such as the High Energy Storage Ring (HESR) planned at 
the GSI laboratory. Quantitative predictions for the charm cross section at such a facility 
are needed in order to assess the prospects for studying charm flavor physics. 

There have been very few previous estimates of the charm production rate by a low-energy 
antiproton beam. Braaten estimated the exclusive cross section for D*°D° in low-energy pp 
collisions by scaling measured cross sections for K*~K + |5j. He obtained a surprisingly 
large cross section of 360 nb at a center-of-mass energy of 5.71 GeV. Titov and Kampfer 
have used a Regge model based on baryon exchange to predict the differential cross sections 
and longitudinal asymmetries for the exclusive production of pairs of charm mesons 0], 
but they did not calculate cross sections integrated over angles. In this paper, we use the 
parton model and perturbative QCD to calculate the inclusive charm cross section from a 
low-energy antiproton beam. 

The parton model and perturbative QCD have been used extensively to calculate cross 
sections for charm quarks, bottom quarks, and top quarks in hadron collisions. In the case of 
charm production, the applications have ranged from fixed-target experiments with proton 
beams with center-of-mass energy y/s below 10 GeV to Fermilab's Tevatron pp collider and 
CERN's Large Hadron Collider with y/s = 2 TeV and 14 TeV, respectively. The parton 
model predictions for top quark production and for the transverse momentum distributions 
for bottom quarks and charm quarks are well tested at the Tevatron. The parton model was 
used to predict the charm cross section in pp collisions at y/s = 200 GeV at Brookhaven's 
Relativistic Heavy Ion Collider (RHIC) [lOj ■ We will use the parton model to make the much 
larger extrapolation down to pp collisions with center-of-mass energies as low as 5.5 GeV. 

A convenient source of flavor-tagged charm mesons in high energy physics experiments is 
the production of the charm mesons D* ± . In the parton model, the inclusive cross section 
for producing D* + in pp collisions can be expressed as 



<t\P(Pa)p(Pb)^D*+ + X] = V / dx x f i/p (x x ) [ dx 2 f. 
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x a[i{x 1 P A )j{x 2 P B ) ^c + X] P C ^ D * + , (1) 

where P c ^d*+ is the fragmentation probability for a charm quark to hadronize into a D* + 
meson. At low energies, the dominant contributions to the sums over the types of partons 
come from = (u,u) and (d, d). The parton distributions fi/p and fj/ p depend on 

a factorization scale [if. The parton cross section a can be calculated as a power series 
expansion in a s (fi r ), with coefficients that depend on /iy and on the renormalization scale 

The lowest energy experiment in which the fragmentation probability P c ^d*+ has been 
measured accurately is e + e~ annihilation at a center-of-mass energy of 10.6 GeV. The Belle 
Collaboration has made precise measurements of the inclusive cross sections for D°, D + , 



Dj, Aj, D , and D* + [151 ] . Using the constraint that the fragmentation probabilities for 
D°, D + , Df, and A+ add up to 1, we can determine the fragmentation probability for D* + 
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at that energy scale to be P c ^d*+ = 24%. We will use this as an estimate of P c ^o*+ in 
low-energy pp collisions. 

To assess whether the extrapolation of the parton model calculation to center-of-mass 
energies as low as 5.5 GeV is plausible, we first consider charm production in e + e~ annihi- 
lation. From the charm meson threshold at y/s = 3.73 GeV up to about 4.6 GeV, the cross 
section varies dramatically due to the charmonium resonances -0(3730), ■0(4040), ■0(4160), 
and ■0(4415). Only above 4.6 GeV is the cross section smooth enough to agree reasonably 
well with perturbative QCD. However quark-hadron duality can be exploited to extend the 
region of applicability of perturbative QCD to significantly lower energies by smearing over 



the energy |14J]. The duality holds between the perturbative charm quark cross section and 
the sum of the cross sections for charm hadrons and for charmonium states below the charm 
threshold. Since the perturbative cross section is smooth without any local minima, the 
smeared hadron cross section should also be smooth. The minimum smearing interval is 
the spacing between radial excitations of the charmonium states, which is about 400 MeV. 
For the perturbative charm quark cross section to be a good approximation to the charm 
hadron cross section alone, the energy must be large enough that the smeared cross section 
receives no significant contributions from the charmonium states below the charm threshold 
and in particular from the -0(3686). This requires the energy at the center of the smearing 
interval to be greater than about 3.9 GeV. 

We now compare inclusive charm production in pp collisions and e + e~ annihilation. One 
difference is in the charmonium resonances that can be produced. The annihilation of e + e~ 
can only produce resonances with quantum numbers 1 . In contrast, pp collisions can 
produce charmonium resonances with any of the possible quantum numbers J PC . Potential 
models predict the spacing between radial excitations in each angular momentum channel 
to be about 400 MeV. Thus smearing in the pp center-of-mass energy by that amount might 
be sufficient to get a smooth cross section above the DD threshold in each J PC channel. 
The sum over J PC will help to further smooth out the cross section. The parton model 
suggests that it may not be necessary to smear over the pp center-of-mass energy, because 
the integration over the momentum fractions x\ and Xi of the colliding partons automatically 
provides some smearing in the energy. To get a smooth cross section, it may be sufficient to 
demand that this integration provides smearing in the center-of-mass energy of the colliding 
partons by at least 400 MeV. 

There is an important difference between e + e~ annihilation and pp collisions in the parton 
processes that create the cc pair and therefore in the hadronization processes that produce 
the charm mesons. The annihilation of e + e~ creates a cc pair and gluons whose overall 
color state is a singlet. In pp collisions, the colliding partons have color connections to 
the remnants of the p and p. The cc pair and the gluons that are created by the collision 
therefore also have color connections to the remnants of the p and p. Part of the energy 
associated with the constituent mass of the light quark in the charm meson may be provided 
by those remnants. Thus it is possible for a parton collision with center-of-mass energy 
below the charm meson threshold to produce charm mesons. 

Based on both the analogy with e + e~ annihilation as well as the differences, we pro- 
pose the following two conditions as plausible criteria for the validity of the parton model 
calculation of charm production: 

A. The integration over parton momentum fractions should provide smearing of at least 
400 MeV in the invariant mass \fl of the colliding partons. 
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FIG. 1: Cross section for inclusive charm production in pp collisions as a function of the center- 
of-mass energy yfs. The solid line is the NLO result for m c = 1.5 GeV and fi r = [if = 2m c . The 
shaded band corresponds to varying m c from 1.3 and 1.7 GeV. The dashed line is the LO result. 



B. The dominant contributions should come from a/J above the DD threshold at 3.73 GeV. 

If the center-of-mass energy is too low for the criteria A and B to be satisfied, the parton 
model calculation can still be used as a phenomenological model for extrapolating the charm 
cross section to low energies. 

To calculate the charm cross section in Eq. ([1]), we use the well-known next-to- leading 
order (NLO) calculations of the parton cross sections for heavy-quark production 0, M, HI- 
The parameters in the charm cross section are the charm quark mass and the renormaliza- 
tion and factorization scales. We do not carry out any resummation of logarithms of the 
transverse momentum, because these logarithms are not large at this low energy. We also 
do not carry out any resummation of the threshold logarithms of s — Am 2 c [la Il7j ] , where 
s is the invariant mass of the colliding partons, although these logarithms can be large. 
The inclusive cross section for pp — > cc + X is calculated using the Monte-Carlo program 
MCFM jllj. The implementation of the cross section at NLO makes use of the matrix ele- 
ments calculated in Ref . 0] • This computation is carried out in the MS scheme with nf = 3 
light flavors of quarks in the initial state. At LO and NLO, we use the parton distribution 
sets CTEQ6L1 and CTEQ6M, respectively 12 ] ■ These sets have been constructed for ri{ = 5 



initial flavors, which leads to a suppression of the gluon density compared to the case rif = 3. 
At NLO accuracy, the correction term that must be added to compensate for the different 
number of flavors is proportional to the gluon-initiated Born-level cross section and depends 
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logarithmically on the ratio /i//m c 13]. Since our factorization scale is comparable to m 
the correction term proves to be small and can be disregarded compared to other theoretical 
errors. 

The parameters m c , fj, r , and fif and their theoretical uncertainties are determined as 
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m c =1 .5 GeV 
m c <^i rf <4m c 
0.5<(a, r /(a, f <2 
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FIG. 2: Cross section for inclusive charm production in pp collisions as a function of the center- 
of-mass energy y^. The solid line is the NLO result for m c = 1.5 GeV and fi r = fif = 2m c . The 
shaded band corresponds to varying fi r and fif from m c to 4m c with ^ < fif/fi r < 2. The dashed 
line is the LO result. 




follows. We choose their central values to be m c = 1.5 GeV and fi r — /if — 2m c . Our 
central value of fi r = fif is twice that used in Ref. M, because this choice decreases the 



sensitivity to the scales at the low energies we consider. The scale 2m c is a natural choice, 
because this is the minimal invariant mass of the virtual gluon in the dominant leading-order 
parton process qq — > cc. We obtain the error band from the charm quark mass by setting 
fi r = fif = 2m c and allowing m c to range from 1.3 to 1.7 GeV. We obtain the error band 
from the scales by setting m c = 1.5 GeV and allowing fi r and fif to range from m c to 4m c 
with ~ < /////i r < 2. 

Our results for the inclusive charm cross section are shown in Figs. [T] and [2] as functions 
of the center-of-mass energy ^fs. Two values of ^/i that are of particular interest are 
4.11 GeV, which corresponds to the maximum antiproton momentum of 8 GeV in Fermilab's 
accumulator ring, and 5.47 GeV, which corresponds to the maximum antiproton momentum 
of 15 GeV in the design of HESR. The central value of the cross section for m c = 1.5 GeV 
and /i r = fif = 2m c at = 5.47 GeV is 12 nb. The error band in Fig. [T]from varying the 
charm quark mass corresponds to a multiplicative factor of about 7.4 2 w 55. The error band 
in Fig. [2] from varying fi r and fif corresponds to a multiplicative factor of about 3.1 2 ~ 9.6. 
The error band from varying m c is significantly wider. 

The reason for the large sensitivity to the charm quark mass is illustrated in Fig. [31 
which shows the distributions at NLO of the invariant mass \fl of the colliding partons for 
\/s = 5.47 GeV and m c = 1.3, 1.5, and 1.7 GeV. The full widths at half maximum of these 
distributions is about 500 MeV, which is large enough to provide the smearing required 
by our criterion A for the validity of the parton model calculation. However the peaks of 
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FIG. 3: Distribution of the invariant mass y§ of the colliding partons at NLO for yfs = 5.47 GeV, 
fi r = fif = 2m c , and m c = 1.3, 1.5, and 1.7 GeV (left, middle, and right curves). The solid vertical 
line marks the position of the DD threshold. 



the distributions are significantly below the DD threshold, which is indicated by the solid 
vertical line in Fig. [3j Thus our criterion B is not satisfied. This casts some doubt on the 
validity of the calculation for y/s as low as 5.47 GeV. 

We can still regard the parton model calculation as a phenomenological model for the 
extrapolation of the charm cross section to low energies. To reduce the model to one with a 
single parameter, we fix the scales at \i r — {if — 2m c and use the charm quark mass m c as 
a phenomenological parameter. One way to determine m c is to require that parton model 
calculations give the measured charm cross sections in low-energy fixed target experiments. 
The most recent data on inclusive charm production from a proton beam on a nuclear target 
with pp center-of mass-energy below 45 GeV are from the LEBC-MPS, ACCMOR, E743, 
E653, E769 and HERA-B Collaborations [HI, [H, H3, HH, H2, HI] - 1 Their beam energies range 



from 200 up to 920 GeV, which corresponds to pp center-of-mass energies from 19.4 up to 
41.6 GeV. The six data points are plotted in Fig. H] as a function of the center-of-mass energy 
yfs. Also shown in Fig. H] is the NLO calculation of the inclusive charm cross section as a 
function of \fs. The solid line is for m c = 1.5 GeV, while the shaded band is the range as 
m c varies from 1.3 to 1.7 GeV. For yfs = 19.4 GeV and m c = 1.5 GeV, the distribution of 



In Ref. 24 1, the various charm meson cross sections published in Refs. 0, llS, 20, 21, 22 1 are converted 
into a common cross section which is denoted by a c s but is actually the inclusive cross section for D 
and D + or, equivalently, the inclusive cross section for D and D~ . To obtain the inclusive charm cross 
section, this must be divided by the sum of the fragmentation probabilities for c — > D° and c — > D + , 
which is approximately 0.81. 
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FIG. 4: Cross section for inclusive charm production in pp collisions as a function of the center- 
of-mass energy y/s. The solid line is the NLO result for m c = 1.5 GeV and fi r = fj,f = 2m c . The 
shaded band corresponds to varying m c from 1.3 and 1.7 GeV. The dashed line is the LO result. 
The data points are the measurements from Refs. 18, 19, 2(], 



21 



22 



23|]. 



the invariant mass \f§ of the colliding partons has a peak near 4.1 GeV and a full width 
at half maximum of greater than 2 GeV. Thus our criteria A and B for the validity of the 
parton model calculation are easily satisfied at this energy. The six data points in Fig. H] 
are more compatible with the smaller values of m c . The best fit to the six data points is 
m c = 1.29 GeV. 

We proceed to use our results to assess the prospects for the study of charm meson mixing 
at HESR. The maximum antiproton momentum is expected to be 15 GeV, which corresponds 
to a center-of-mass energy of 5.47 GeV. The NLO parton model with m c = 1.29 GeV gives 
a cross section of 89 nb at this energy. The design luminosity of HESR is 2 x 10 32 /(cm 2 s). 
Using the NLO cross section, we find that 5.6 x 10 8 charm events could be produced in a year 
of dedicated running. Using our estimate P c ^d*+ = 24% for the fragmentation probability, 
we obtain 2.7 x 10 8 D* + or D*~ events. The tagging of the initial flavor of the charm meson 
as D° or D° comes from the decay into D°tt + or D°tt~, which has a branching fraction 
of 68%. The simplest decay mode with which to observe mixing is the doubly-Cabbibo 
suppressed decay mode D° —>■ K + ir~ , which has a branching fraction 1.5 x 10~ 4 . Thus our 
estimate of the cross section corresponds to 27, 000 D° — > K + tt~ events. To put this into 
context, the evidence for charm meson mixing at the Tevatron from the CDF Collaboration 
used 13,000 D° — > K + n~ events [4J. We conclude that the charm cross section at HESR 
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could be large enough to study charm meson mixing. 

There is also a proposal to use the antiproton source for Fermilab's Tevatron for low- 



energy pp collisions [25] . The maximum antiproton momentum in the accumulator ring at 
Fermilab is 8 GeV, which corresponds to a center-of-mass energy of only 4.11 GeV. The 
curves in Figs. [T] and [2] have been extended down to the DD threshold at 3.73 GeV, but 
this is too low for the model to be plausible, because the cross section must vanish at 
this threshold. A center-of-mass energy of 4.11 GeV may also be beyond the domain of 
plausibility of the model, because it is only a little above the D*D* threshold. If we ignore 
this problem, the NLO parton model with m c = 1.29 GeV gives an inclusive charm cross 
section of 4.1 nb. 

In summary, we have used the parton model and NLO perturbative QCD to calculate 
the cross section for inclusive charm production by a low-energy antiproton beam. At y/s = 
5.47 GeV, which corresponds to a p beam with momentum 15 GeV, the cross section can be 
increased or decreased by an order of magnitude by varying m c , /i r , and /if within reasonable 
ranges. Since the fundamental justification for the calculation is questionable at such low 
energies, we have proposed the parton model calculation as a phenomenological model for 
extrapolating the charm cross section. Treating m c as a phenomenological parameter, we 
determined its value by fitting measured charm cross sections from fixed-target experiments 
with low-energy proton beams. The resulting estimate of the charm cross section at \/s = 
5.47 GeV is 89 nb. This cross section is large enough to allow the study of charm meson 
mixing in experiments with a 15 GeV antiproton beam. 
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